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ABSTRACT 

We present new 7 mm and archive 1.3 cm high angular resolution observations of the HL/XZ Tau region 
made with the VLA. At 7 mm, the emission from HL Tau seems to be arising in a clumpy disk with radius of 
order 25 AU. The 1.3 cm emission from XZ Tau shows the emission from a binary system with 0."3 (42 AU) 
separation, known from previous optical/IR observations. However, at 7 mm, the southern radio component 
resolves into a binary with 0."09 (13 AU) separation, suggesting that XZ Tau is actually a triple star system. 
We suggest that the remarkable ejection of gas from the XZ Tau system observed with the HST may be related 
to a periastron passage of this newly discovered close binary system. 

Subject headings: ISM: individual (XZ Tau, HL Tau) — ISM: jets and outflows — radio continuum: ISM — 
stars: formation 
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1. INTRODUCTION 

The region encompassing the stars HL/XZ Tau, in the 
northeastern part of the LI 551 dark cloud, has been the sub- 
ject of many studies over the years. This region lies at a dis- 
tance of 140 pc (e.g., Kenyon et al. 1994; Torres et al. 2009) 
and it is particularly rich in HH jets, being one of the regions 
where this phenomenon was first identified (Mundt & Fried 
1983). 

HL Tau is one of the most intensively studied T Tauri stars. 
Since the first proposal that this star is associated with a nearly 
edge-on circumstellar disk (Cohen 1983), numerous studies 
have been carried out in order to image the proposed disk 
(e.g., Sargent & Beckwith 1991 ; Wilner et al. 1996; Looney et 
al. 2000). This star has been proposed as the exciting source 
of a molecular outflow (e.g., Torrefies et al. 1987; Monin et 
al. 1996). HL Tau is also the source of a collimated optical 
jet-counterjet system that has been extensively studied (e.g., 
Mundt et al. 1990; Rodriguez et al. 1994; Anglada et al. 
2007). Very recently, Greaves et al. (2008) presented high 
angular resolution maps at 1.3 cm of the jet-disk system in 
HL Tau. These authors reported the detection of a weak 1.3 
cm source which they propose is a 14 Mj proto-planet orbit- 
ing at a radius of ~65 AU around HL Tau. 

XZ Tau, located ~25" to the east of HL Tau, is a close 
binary system composed of a T Tauri star and a cool compan- 
ion separated by 0."3 (Haas et al. 1990). XZ Tau is also the 
source of an optical outflow, as revealed, e.g., by the studies 
of Mundt et al. (1990). The spectacular sequence of Hubble 
Space Telescope (HST) images of Krist et al. (1999) shows 
evidence of the expansion of nebular emission, moving away 
from XZ Tau with a velocity of ~70 km s -1 . This bubble of 
nebular emission is different from the collimated jets seen in 
other young stars, for example HL Tau. Recent new HST ob- 
servations (Krist et al. 2008) revealed a succession of bubbles 
and a fainter counterbubble. In addition, they reveal that both 
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components of the XZ Tau binary are also driving collimated 
jets. Krist et al. (2008) proposed that the bubble is the result 
of a large velocity pulse in the collimated jet driven by the 
southern component of the XZ Tau system. 

In this Letter, we present the results of sensitive high- 
angular resolution Very Large Array (VLA) observations at 
1.3 cm and 7 mm of HL/XZ Tau. These observations allowed 
us to map the dust emission from the disk of HL Tau at scales 
of ~7 AU and reveal that XZ Tau is actually a triple-star sys- 
tem with the southern optical component resolved into a close 
binary with a separation of only 13 AU. 

2. OBSERVATIONS 

Observations at 7 mm continuum were made using the VLA 
of the National Radio Astronomy Observatory (NRAO) 4 . The 
data were taken in A, B and D configurations (VLA Project 
Codes: AC0763, AC0850 and AC0816). Total on-source in- 
tegration times were ~3 h at A Configuration (observed in 
three runs during 2004 October 15, 17 and 26), ~1 h at B 
Configuration (2006 September 19), and ~1.5 h at D Con- 
figuration (observed in two runs during 2006 January 20 and 
22). The fast- switching technique was used for all the ob- 
servations. Phase and flux calibrators were 0431+175 and 
3C286, respectively. The same phase center was used for all 
the observations. 

We also used archive continuum data at 1.3 cm taken with 
the VLA in its A Configuration (VLA Project Code: AG071 1) 
with a total on-source integration time of ~12 h (observed in 
two runs during 2006 March 12 and April 1). Phase and flux 
calibrators were 0431+206 and 3C286, respectively. The Pie 
Town antenna was used originally in these observations, but 
we did not include it in our analysis. 

Data editing and calibration were carried out using the 
Astronomical Image Processing System (AIPS) package of 
NRAO, following the standard high-frequency VLA proce- 
dures. Primary beam corrections were applied using the task 
PBCORofAIPS. 

Maps at 7 mm were made by concatenating the data from all 
the configurations (A, B and D). This allowed us to improve 
by ~20% the signal-to-noise ratio of the A configuration data 

4 The NRAO is a facility of the National Science Foundation operated 
under cooperative agreement by Associated Universities, Inc. 
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Fig. 1 . — (a) Superposition of the 7 mm map (contours) over the 1.3 cm continuum map (color scale) of HL Tau. Contours at 7 mm are -3, 3, 4, 5, 6, 8, 10, 
12, 14, 18, 22, 26 and 32 times the rms of the map, 50 fjjy. The dashed lines mark the directions of the optical jet (P. A =50°) and the disk (P. A =140°). The half 
power contour of the 1.3 cm (filled ellipse) and 7 mm (empty ellipse) beams are shown in the bottom left corner, (b) A 7 mm intensity profile along the major 
axis of the disk (P.A.=140°), obtained by averaging the emission at ±0"05 from the axis marked in (a). The intensity error is indicated by the thickness of the 
profile. The disk truncates at a radius of ~12 AU to the NW, and shows a gap at a radius of ~10 AU to the SE. We speculate that these features could be related 
to the formation of a proto-planet (see text). 



without significant detriment of the angular resolution. The 
7 mm map of HL Tau shown here was made using data with 
baselines >150 kA in order to resolve out extended emission 
that surrounds HL Tau. The nature of this extended emis- 
sion is not fully understood, but probably originates from dif- 
fuse dust in the region. The 1.3 cm map of HL Tau shown in 
Greaves et al. (2008) was obtained from the same data, and, 
as expected, is very similar to our 1.3 cm map. 

The rms noises of the natural weighting maps were 13 /iJy 
at 1.3 cm, and 50 /iJy at 7 mm. Synthesized beams are ~0."10 
at 1.3 cm and ~0."05 at 7 mm. Bandwidth smearing at the 
position of XZ Tau degrades the resolution by ~5% in the 
radial direction, which does not affect our analysis. 

This region has a known proper motion of 0. // 025±0. // 0035 
yr -1 (Rodriguez et al. 2003), so we shifted the 2006 1.3 cm 
data by 0."04 to align them with the 2004 7 mm data (with 
uncertainties of 0."005 in proper motions and O'.'Ol in absolute 
astrometry at each epoch). 

3. RESULTS AND DISCUSSION 
3.1. HL Tau 

In Figured we show a superposition of the 7 mm emission 
of HL Tau (contours) over the 1.3 cm map (colors). At 1.3 
cm, HL Tau shows a quadrupolar morphology consisting of 
two nearly perpendicular structures at position angles (P.A.s) 
of ~50° and ~140°. As previously noted by Greaves et al. 
(2008), the alignment of the 1.3 cm emission at P.A.~50° 
with the optical jet of HL Tau (e.g. Anglada et al. 2007), 
suggests that the 1.3 cm emission at this P. A. is tracing free- 
free emission from the shock-ionized gas at the base of the 
flow. In contrast, the 1.3 cm emission at P.A.~140° is most 
probably tracing thermal dust emission from a circumstellar 
disk perpendicular to the outflow. This is more evident in our 



7 mm map (see Fig. [T^), where the emission shows an elon- 
gated morphology with a PA. of ~140° (perpendicular to the 
jet) and a size of ~50 AU. The positions of the 1.3 cm and 7 
mm maxima are nearly coincident and define the center of the 
radio jet. 

We also detected two weak compact 1.3 cm features, one 
~0."4 to the NW (at a level of 3-cr) and the other ~0."4 to 
the S (at a level of 4-cr) of the 1.3 cm emission peak (see Fig. 
Ek). Greaves et al. (2008) identified the 1.3 cm NW feature 
with a 1.4 mm source detected by Welch et al. (2004), ob- 
taining a spectral index of 2.5 in this wavelength range. They 
proposed that the 1.3 cm NW feature is dust emission arising 
from a proto-planet with a mass of 14 Mj and orbiting around 
HL Tau at a radius of ~65 AU. The other 1.3 cm feature to 
the S, although somewhat stronger than the NW feature, was 
not discussed by these authors. As can be seen in Figure Hk, 
we have not detected 7 mm emission at the position of the 
proto-planet candidate, with a 3-cr upper limit level of ~0.15 
mJy. This upper limit and the 1.3 cm flux density of Greaves 
et al. (2008) implies a spectral index <1, suggesting free-free 
rather than dust emission. On the other hand, the spectral in- 
dex of 2.5 between 1.3 cm and 1.4 mm obtained by Greaves 
et al. (2008) would imply a 7 mm flux density of 0.5 mJy, 
well above our 3-cr upper limit. Additional, high angular res- 
olution observations in the mm range are necessary to clarify 
the true nature of this source. 

In Figure[l}) we show an intensity profile of the 7 mm emis- 
sion along the major axis of the disk. As can be seen in Fig- 
ures andHJ), the 7 mm emission is knotty and asymmetric. 
The proposed disk extends ~25 AU to the SE, and shows a 
gap at a radius of ~10 AU and the emission to the NW ap- 
pears truncated at a radius of ~12-15 AU. The gap and the 
truncation are probably related to a decrease of the density in 
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Fig. 2. — Radio continuum maps of XZ Tau at 3.6 cm (a), 1.3 cm (b) and 7 mm (c). The contours are -3, 3, 4, 5, 6, 8, 10, 12, 14, 18, 22 and 26 times the rms 
of each map, 7 /iJy (3.6 cm), 13 /iJy (1.3 cm) and 50 /iJy (7 mm). At 3.6 cm a single radio source is detected. The 1.3 cm map shows two continuum sources, 
labeled as N and S. The source S is marginally resolved, showing a slight elongation to the south. Crosses mark the positions of the two components SI and S2 
detected at 7 mm. At 7 mm, the source N is detected as a weak source, and the 1.3 cm source S, is resolved in two components, SI and S2. 



TABLE 1 

Components of the XZ Tau triple system 



Radio Optical Position (J2000) a Flux Density (mJy) a Spectral Emission 

Component Component RK DEC Error 1.3 cm 7 mm Index Mechanism 

N XZ Tau B 04 31 40.0751 18 13 57.079 77 006 0.14 ±0.03 0.3 ± 0.2 1.1 ± 1.0 Free-free 

51 XZ Tau A 04 31 40.0891 18 13 56.853 0"005 0.33 ± 0.03 b 1.2 ± 0.2 1.9 ± 0.3 Free-Free + Dust 

52 XZ Tau C c 04 31 40.0899 18 13 56.762 0"008 0.08 ± 0.03 b 0.8 ± 0.2 3.4 ± 0.7 Dust 



Notes.- Units of right ascension are hours, minutes, and seconds, and units of declination are degrees, arcminutes, and arcseconds. 

a Derived from Gaussian fits to the 7 mm natural weighting map (Epoch 2004.8). Errors in positions correspond to the relative positions. The error in the absolute position is estimated 
to be r^o'/oi^ Flux density derived by fitting the marginally resolved 1.3 cm source S to two unresolved Gaussian sources. This component is probably a deeply embedded object 
undetectable at optical/IR wavelengths. 



the disk. Thus, since the disk is observed nearly edge-on (e.g. 
Cohen 1983), this gap could be associated with an annulus of 
low density in the disk with a radius of ~1 0-1 5 AU, similar to 
the radius of the orbit of Saturn, or twice the orbit of Jupiter. 
Such structures are expected to form as a consequence of the 
accretion of a planetary object forming in the disk (see, e.g. 
Bryden et al. 1999). Therefore, the gap and the truncation 
of the disk detected in our 7 mm map could be related to the 
formation of a proto-planet. Alternatively, this gap could be 
due to tidal truncation in a binary system (e.g., Pichardo et al. 
2005) as in the case of HH30 (Anglada et al. 2007, Guilloteau 
et al. 2008). However, with the present data, these sugges- 
tions are very speculative, and should be confirmed by subse- 
quent observations and/or modelling of the disk of HL Tau. 

3.2. XZTau 

In Figure^ we show the 3.6 cm continuum map (angu- 
lar resolution ~0."3) of XZ Tau, obtained by Rodriguez et 
al. (1994), where a single radio continuum source was de- 
tected. Figures [2j) and Et show our 1.3 cm and 7 mm maps, 
respectively. At 1.3 cm (angular resolution ~0."1), XZ Tau 
is resolved in two radio continuum sources (labeled N and S) 
separated by ~0."3 with a P.A. of ~143°. The 7 mm map (an- 
gular resolution ~0."05), reveals that the 1.3 cm component S 
is actually a double radio source (that we label as components 
SI and S2) with a separation of 0."09 and a P.A. of -6°. 

The size of the synthesized beam of the 1.3 cm map (~0."1), 
slightly larger than the separation between the S 1 and S2 com- 
ponents at 7 mm (~0."09), does not allow to completely sep- 



arate the emission of these components at 1.3 cm. However, 
the 1.3 cm source S is marginally resolved in our map, show- 
ing a slight elongation to the south. In Figure [2J3 we show 
the positions of the components S 1 and S2 derived from the 
7 mm map. As can be seen, most of the 1.3 cm emission 
seems to arise from component SI, since this is the source 
nearest to the peak of the 1.3 cm source S. The 1.3 cm source 
S can be fitted by two unresolved Gaussian sources separated 
by ~0."08, very similar to the separation between sources SI 
and S2 in the 7 mm map. 

In Table [T] we give the position (derived from the 7 mm 
map), flux density at 1.3 cm and 7 mm, and the spectral index 
for each component. For source N, we obtained a spectral in- 
dex of 1.1 ± 1.0 consistent with free-free emission from this 
source. For source S as a whole, we obtained a "combined" 
spectral index of 2.5 ± 0.3 (obtained from the total flux den- 
sity of the 1.3 cm source S, and the sum of the 7 mm flux 
densities of SI and S2). This steep spectral index suggests 
a combination of free-free plus thermal dust emission. If we 
use the flux densities obtained from a double Gaussian fit to 
the 1.3 cm source S, for source SI we obtain a spectral index 
of 1.9 ± 0.3 which suggests a combination of free-free plus 
thermal dust emission. In contrast, for S2 we obtain a spectral 
index of 3.4 ± 0.7 which suggests that the emission of this 
source is dominated by thermal dust emission. 

As mentioned above, XZ Tau is known to be an optical/IR 
binary system with a separation of ~0."3. Following the nota- 
tion of Krist et al. (2008), we will refer to the optical compo- 
nents of this system as XZ Tau A (the southern optical compo- 
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Fig. 3.— (a) Plot of the P. A. between the optical/IR components XZ Tau A and XZ Tau B vs. time (Haas et al. 1990; Krist et al 1999, 2008). The dashed 
line is a least squares fit to these data, and corresponds to an orbital angular velocity of -0.9 ±0.1 degrees yr -1 . Three different P.A.s measured from our 1.3 cm 
and 7 mm maps are marked with diamonds. Only the P. A. between the 7 mm sources N and SI seems to be compatible with the angular orbital velocity derived 
from the optical/IR data. This suggests that the source SI is the radio counterpart of the optical component XZ Tau A. (b) Plot of the separation between the 
optical/IR components XZ Tau A and XZ Tau B vs time. The measurements are compatible with a constant separation of 0"30±0"01. (c) and (d) A schematic 
representation of the XZ Tau triple system. The positions of XZ Tau B relative to XZ Tau A, derived from the P. A. and separation measurements of Haas et al. 
(1990), Krist et al (1999, 2008) and our data, are shown. The current data are compatible with a circular face-on orbit of XZ Tau B relative to XZ Tau A. In 
contrast, with only a single detection of the third component, XZ Tau C, we have no information about its orbital motion. Two hypothetical orbits (circular and 
elliptical) for XZ Tau C are shown. 



nent) and XZ Tau B (the northern optical component). Krist 
et al. (1999, 2008) observed XZ Tau with the HST in several 
epochs, and detected a change in the P.A. of the binary due to 
orbital motions. According to this motion (~ -0.9 deg yr -1 ) 
at the epoch of our radio observations, component XZ Tau B 
must be still located north of component XZ Tau A. There- 
fore, we identify the radio source N with the optical compo- 
nent XZ Tau B. 

The identification of the two radio sources at the south (S 1 
and S2) with the other optical component (XZ Tau A) is not 
straightforward. One possibility is that both radio sources 
have optical emission that could not be resolved by the HST 
observations of Krist et al. (2008). However, another pos- 
sibility is that only one of the radio sources is the counter- 
part of the optical component. In order to investigate these 
possibilities, we have plotted in Figures [3^ and [3J) the P.A. 
and separation of the optical binary measured in each epoch 



together with the P.A.s and separations between the differ- 
ent radio components measured from our 1.3 cm and 7 mm 
maps. As can be seen, the P.A. and separation between the 
radio sources N and S 1 seems to be the only one compatible 
with the previous optical/IR data. This suggests radio source 

5 1 as the sole counterpart to the optical component XZ Tau A; 

52 is optically undetected. This is similar to the case of the 
close binary SVS 13, where both components are detected at 
cm wavelengths, but only one object is detected in the optical 
regime (Anglada et al. 2004). As we have discussed above, 
the radio emission of component S2 seems to be dominated 
by thermal dust emission. We believe that this source traces a 
star since if it not were gravitationally bound, it would quickly 
dissipate. The size of Sb is smaller than 7 AU, that expanding 
at only 1 km s" 1 , will dissipate in only 30 years. Then, we 
conclude that this dust must be gravitationally bound to a star 
and it is most possibly an envelope or a disk. This suggests 
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that source S2 is a more deeply embedded object so it is not 
surprising that it is optically very faint. Therefore, we propose 
that the radio source S 2 is an optically invisible third compo- 
nent of the XZ Tau system. For consistency with the previous 
nomenclature, we call this component XZ Tau C. 

In summary, the results derived from our 1.3 cm and 7 mm 
maps suggest that XZ Tau, previously known as a binary sys- 
tem, is actually a triple star-system (see Figs. [3J: and Ol). 
The previously known optical/IR components, XZ Tau A and 
XZ Tau B, have 1.3 cm and 7 mm counterparts and are sepa- 
rated by 42 AU. A third component, XZ Tau C has been de- 
tected for the first time in our 1.3 cm and 7 mm maps, and is 
found to be a more deeply embedded object, probably unde- 
tectable at optical/IR wavelengths, that forms a close binary 
system with XZ Tau A with a separation of only 13 AU. 

Each optical component of the XZ Tau system (A and B) is 
known to drive a collimated jet (Krist et al 2008). It has been 
proposed that the jet of XZ Tau A underwent a large velocity 
pulse approximately in 1980 (Krist et al. 2008) which, after 
shocking and heating the gas near the star, would result in the 
70 km s -1 expanding bubble detected in the series of HST im- 
ages of Krist et al. (1999, 2008). These authors explored the 
possibility that this ejection was triggered by the passage of 
component XZ Tau B very near component XZ Tau A. How- 
ever, even assuming a very eccentric orbit (e^0.9) they de- 
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rived that the last approach of these components would have 
taken place in 1955, and thus, they discarded this possibility. 

The discovery of a new companion, closer to XZ Tau A, 
suggests that it is once more feasible to consider the passage 
of this companion as the ejection trigger. If we assume a face- 
on, nearly circular (e^O.l) orbit and a total mass of 1 M 
(as suggested by the results of Krist et al. 1999 and White 
& Ghez 2001), we derive an orbital period of ~40 yr. If at 
the time of our 7 mm observations (2004.8) the XZ Tau A- 
C system wass near apoastron, then the last periastron pas- 
sage must have taken place one half-period before, i.e. nearly 
1980, the date estimated for the ejection. Further observa- 
tions are needed to produce a more accurate ephemeris for 
the XZ Tau system. The present results predict a new perias- 
tron passage around 2020, and then a new major ejection from 
XZ Tau A at this epoch. 
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